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Abstract. We propose a geometric perspective to describe the motion 
of self-propelled particles moving at constant speed in d dimensions. We 
exploit the fact that the vector that conveys the direction of motion of 
the particle performs a random walk on a {d— l)-dimensional manifold. 
We show that the particle performs isotropic diffusion in d-dimensions 
if the manifold corresponds to a hypersphere. In contrast, we find that 
the self-propelled particle exhibits anisotropic diffusion if this manifold 
corresponds to a deformed hypersphere (e.g. an ellipsoid). This simple 
approach provides an unihed framework to deal with isotropic as well 
as anisotropic diffusion of particles moving at constant speed in any 
dimension. 


1 Introduction 

We find examples of self-propelled entities in a remarkably large variety of chemi¬ 
cal [1-6], physical [7-9] and biological [10] systems. The non-equilibrium nature of 
these active systems has recently been investigated through the development of what 
has been called active soft matter [10,11]. In active systems, such as self-propelled 
particle (SPP) systems, the constant conversion of energy into work - used by the 
particles to self-propel in a dissipative medium - drives the system out of thermo¬ 
dynamic equilibrium leading to remarkable physical properties for both, interacting 
as well as non-interacting SPP systems. In interacting SPP systems, large-scale col¬ 
lective motion patterns [10-12] and non-equilibrium clustering [13-15] are observed 
in the presence of a velocity-alignment interaction for homogeneous media with peri¬ 
odic boundary conditions. Moreover, the presence of fluctuations in both, the moving 
direction of the particle and its speed - typically related to fluctuations of the self- 
propelling engine - leads to bistability of macroscopic order and disorder [16]. For 
non-interacting active particles, these non-thermal fluctuations result in complex non¬ 
equilibrium transients in the mean squared displacement [17] and anomalous (non- 
Maxwellian) velocity distributions [6,18]. Interestingly, the physics of active systems 
is remarkably different when other boundary conditions are used [19] - the lack of 
momentum conservation in active systems induces non-classical particle-wall interac¬ 
tions which allow, for instance, the rectification of particle motion [20-24]. Besides, 
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non-interacting SPPs can exhibit spontaneous particle trapping and subdiffusion in 
heterogeneous media [25,26]. 

In this work, we focus on the simplest class of SPP systems that we can think of: 
non-interacting SPPs moving at constant speed. Despite their simplicity, these SPP 
models find applications in artificial active particles such as vibration-driven rods and 
disks [7-9], light-driven [27-30] and chemically-driven [1-6] SPPs at low density as 
well as diluted bacterial systems [31]. Here, we contribute to the theoretical descrip¬ 
tion developed in previous works [6,17,18,32] by proposing a geometric perspective 
onto self-propelled motion. A geometric formulation of the model has several advan¬ 
tages. First, geometric insights yield the possibility of understanding observations 
intuitively and allow to generalize two dimensional SPP models to any spatial dimen¬ 
sion immediately. Furthermore, we show that the same theoretical framework enables 
us to describe SPP motion in anisotropic environments as occurs, for instance, in 
experiments with eukaryotic cells moving on a pre-patterned surfaces [33,34]. More¬ 
over, the geometric formulation simplifies the mathematical modeling itself: we de¬ 
rive stochastic equations of motion from a geometric principle which is reflected by a 
Fokker-Planck equation defining an unique stochastic process. Accordingly, ambigui¬ 
ties related to the interpretation of Langevin equations [35] do not occur. In short, we 
propose a simple unified framework to model isotropic as well as anisotropic diffusion 
of SPPs in any dimension. 

This paper is organized as follows. In section 2, we present our geometric approach 
to model self-propelled motion and discuss technical details in section 3. Based on 
this general framework, we review isotropic self-propelled motion in section 4 focusing 
on two spatial dimensions and subsequently generalizing to arbitrary dimensions. In 
section 5, we use the geometric approach to study SPPs in anisotropic environments. 
We summarize and discuss our results in section 6. 


2 A geometric view on self-propelled motion 

We consider non-interacting SPPs moving at constant speed vq. The moving direction 
is denoted by a vector e(t), called director for short in the following. The dynamics 
of an active particle is then given by the following differential equation: 

= v(t) = uoe(t), e{t) =e{t)/\e{t)\. (1) 

The spatial dynamics is solved by 

r(t) = r(to) + voj^ dt' (2) 

In this sense, the spatial dynamics r(t) of the SPP is subordinated to the stochastic 
dynamics of the director e(t). Once the stochastic properties of the director e(t) are 
known, the motion of a SPP in space can be characterized using Eq. (2). For example, 
the correlation function of the director and the mean squared displacement of a SPP 
are linked by the Taylor-Kubo formula [36-38]: 

(k(t) - r(to)l^) = (l^r(f)l^) =Vo [ dt' f dt” (e(t') • e{t”)). (3) 

Jto Jto 

Having understood that the motion of an active particle moving at constant speed is 
prescribed by the dynamics of the director, we focus on the latter henceforward. We 
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Fig. 1. (a) Visualization of the trajectory r{t) of a SPP in three dimensions and (b) the 
corresponding diffusive motion of the moving direction e(t) on a sphere S^. Notice that the 
SPP moves in d dimensions - as shown in panel (a) - and the moving direction e(t) lives on 
a surface of dimension d — 1, panel (b). 


will show that a geometric view onto the director dynamics simplifies the modeling 
of self-propelled motion in arbitrary dimensions, both in isotropic and anisotropic 
environments. 

Let us briefly recall how self-propelled motion is modeled in two dimensions [32,39]. 
Usually, the director is parametrized by an angle (p{t) as follows 

where the polar angle obeys the stochastic differential equation 

( 5 ) 

Hence, the trajectory described by the director e(t) is equivalent to the trajectory 
of a Brownian particle moving on a circle. Accordingly, the two dimensional motion 
of a SPP is subordinated to the ordinary Brownian motion of the director e(t) on 
a manifold that in this case is simply a circle. This geometric interpretation of Eq. 
(4)-(5) enables us to generalize the model to arbitrary dimensions as follows. The 
motion of a SPP in three dimensions is controlled by the director dynamics e(t) which 
describes diffusive motion on a sphere as illustrated in Fig. 1. Thus, the director 
describes Brownian motion on a {d — l)-sphere in the case of d-dimensional self- 
propelled motion. In other words, the description of isotropic self-propelled objects is 
directly related to random walks on hyperspheres [40] where Eq. (2) links these two 
processes. 

If the trajectory e(f) corresponds to Brownian motion on a perfectly isotropic 
sphere, the resulting trajectory r(t) of the SPP is isotropic in space as well. Ac¬ 
cordingly, the SPP does not possess a preferential direction of motion. Following the 
geometric picture drawn above, the description of anisotropic self-propelled motion is 
straightforward: the director e(t) describes Brownian motion on a compact manifold 
different from a hypersphere. Such a manifold could be, for instance, an ellipsoid. 

In summary, the SPP motion is isotropic if the director performs Brownian motion 
on a hypersphere. In contrast, SPP motion is (typically) anisotropic if this manifold 
corresponds to the surface of a “deformed” hypersphere (the details on the properties 
of the manifold are given in section 3). 
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Thus, to close our simple SPP model we have to (i) specify the surface on which 
the director moves and (ii) provide the dynamics for the director on this surface. 
We argued above that the motion of the director on the manifold is purely diffusive. 
Hence, the dynamics of the probability density for the director pointing in a certain 
direction, P(e, t), is determined by the diffusion equation 


dP{e, t) 
Wt 


DAeP{e,t). 


( 6 ) 


The operator denotes the Laplace-Beltrami operator on the particular manifold. 
The parameter D measures the strength of the fluctuations acting on the director 
which defines the diffusion coefficient on the manifold but not the diffusion coefficient 
of the SPP. Notice that according to (6), D determines the characteristic time scale 
of relaxation but, by definition, does not affect stationary states. Depending on the 
particular system under consideration, the noise strength D may depend on addi¬ 
tional parameters, for example the speed of the particle describing speed dependent 
persistence [18,32]. In this study, we will treat it as an independent parameter. 


The spatial diffusion coefficient 

Let us suppose that the underlying stochastic process described by the director is sta¬ 
tionary [35] in the sense that observables do not depend on the instant of time at which 
they have been measured. Under this assumption, the director correlation function 
does depend on time differences only and we may write C{At) = (e(P) • e(P -|- At)). 
We define the correlation time Tc as integral of the correlation function^: 

dAtC{At). (7) 

Using this definition and assuming (e(t —)■ oo)) = 0, one can rewrite integral (3) in 
order to obtain the long time limit 

~ 2Tcvlt = 2d'D^t. 

The last expression in the equation above is the defining relation of the spatial dif¬ 
fusion coefficient Px which is related to the correlation time and the spatial dimen¬ 
sionality via 



V 


X 


TcVl 

d 


( 8 ) 


The preceding discussion illustrates that the stochastic motion on curved space de¬ 
scribed by the director as well as the properties of SPPs are inextricably linked. 


3 Langevin dynamics of the director 

According to the model explained in the previous section, the dynamics of the director 
corresponds to Brownian motion on a compact manifold [40,41]. We assume that this 

^ This definition is valid only if the correlation function does not oscillate. Moreover, 
we implicitly assume that the integral over the correlation function converges excluding 
anomalous diffusion processes in this context. 
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compact manifold can be embedded in d dimensions and parametrize the points of 
this space by Cartesian coordinates: e = (ci, 62 ,e^)- In subsequent sections, we 
will associate the direction of motion of the SPP with a point on the manifold via 
e{t) = e(t)/|e(t)|, cf. Eq. (1). In this section, we derive the stochastic differential 
equations that govern the stochastic motion of the director e on the manifold whose 
dynamics is defined by the Fokker-Planck equation (6). 

At a more mathematical level, the comment above means that the director is a 
d-dimensional vector e G while its motion is confined to a surface of a geometrical 
object (a compact manifold) which is embedded in d dimensions. Let the points on the 
manifold be parametrized by n < d generalized coordinates where /x = 1, 2,..., n. 
For concreteness, imagine the coordinates to be the angles which parametrize 
the points on an unit sphere. The transformation of coordinates from the Cartesian 
laboratory reference frame to generalized coordinates reads 

ei = (9) 


where t = 1,2,..., d. We define the Jacobian matrix J of this transformation as 


_ dcr 

d^i,' 


( 10 ) 


The central objects for the description of the diffusion on a manifold are the metric 
tensor g whose elements read 




E 

i=l 




( 11 ) 


as well as the inverse metric tensor A = g~^. 

In (9), we expressed the director by generalized coordinates. The Laplacian oper¬ 
ator in Eq. (6) must be expressed by these coordinates as well [42] 


I It ‘It 

E E (vidi ) 


( 12 ) 


where 9^ is a shorthand for the derivative with respect to and we denote the 
determinant of the metric tensor by jgj. The probability density P{e,t) becomes a 
function of the generalized variables and it is normalized according to the condition 

J dpi J dp2 ■■■ j dpn \/M t) = 1- (13) 

We introduce the probability density p({(p}, t) to find the director pointing in a certain 
direction by absorbing the measure factor in (13) as 


P{W}A) = V\9\P{e,t)- 


(14) 


The dynamics of p({(/?}, t) follows from its definition and the general form of the 
Laplacian operator on the manifold. We reorganize terms in the resulting equation 
such that it takes the standard form of a Fokker-Planck equation: 


dp{W},t) 

dt 


= py^d^d^ - Py^d^ p (duAf^i, + A^^di,logy/\g\"j . (15) 



6 


Will be inserted by the editor 


From this Fokker-Planck equation which is written in the so called Ito form [35] we 
can immediately read off the corresponding Langevin dynamics for the generalized 
variables. The first term encodes the strength of the fluctuations (diffusion term), 
whereas the second term describes a deterministic drift. In order to write down the 
Langevin dynamics for the generalized coordinates, we introduce the generalized noise 
amplitudes cr, a (n x n)-matrix having the property yl = cr • cr^. Hence, the stochastic 
differential equation for the generalized coordinates in Ito (I) interpretation reads 

(I) +Af,^,^^log^/\g\y (16) 

CX V 

The corresponding differential equation in Stratonovich (S) interpretation is obtained 
by applying the usual rules of stochastic calculus [35]: 

(S) 

^ flCt (^) ^ ^ ^ I E fiOL “ 1 “ ^y^og^/\g\ 

oc \ ct 

(17) 

The random processes (t) denote independent Gaussian random processes with zero 
mean and temporal ^-correlations. Since the dynamics involves multiplicative noise - 
the noise amplitude depends on the state of the system - a word on the interpretation 
of the Langevin equation is in order. Both Langevin equations, (16) and (17), describe 
the same physics since they were derived from the same Fokker-Planck equation. In 
other words, one can choose the interpretation depending on which equation is easier 
to handle analytically or more convenient to use in a numerical experiment. 

The equations (16) and (17) describe diffusion on a compact manifold which is 
embedded in d-dimensional space. They are simplified considerably if the manifold is 
parametrized using an orthogonal basis. In this particular case, the metric tensor is 
diagonal and we find 


and thus 


(I) 

(S) 



0^1, 

(18) 

+ 9^^ log ^/\g\ ), 

(19a) 


(19b) 


In (19), both drift and diffusion terms are due to geometric properties of the 
manifold and proportional to the strength of the fluctuations D. This last point, 
i.e. that drift and diffusion terms are proportional to D, distinguishes our approach 
from the motion of particles subjected to an external forcing. In contrast to motion 
induced by an external force, the drift and diffusion terms are not independent, in 
the sense that by taking the limit Z? —> 0, both terms vanish, while we expect the 
drift term to survive in the presence of an external forcing. In summary, the strength 
of an external field or the response of a particle to an external stimulus requires some 
additional parameter (force strength). For anisotropic environments, we will make 
use of a parameter characterizing the anisotropy of the environment for instance, but 
this will not be equivalent to an external force. 
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In general, the Langevin equations (16) and (17) determine the temporal evolu¬ 
tion of the generalized coordinates given any arbitrary metric tensor. However, we 
want to interpret the vector e as the direction of motion of a self-propelled parti¬ 
cle, cf. Eq. (1) and Eq. (9). Therefore, the embedding of the manifold is not arbi¬ 
trary. In principle, the vector e could be shifted by a constant vector Gq (indepen¬ 
dent of the generalized coordinates) without changing the metric tensor or the equa¬ 
tions of motion. However, the velocity of the SPP would shift according to Eq. (1): 
V = upe —)■ Up (e -b eo)/|e -b eo|. The constant shift would imply ballistic motion into 
one particular direction determined by ep. We fix this constant by requiring that 

Jdifi Jd(p 2 ■ ■ ■ Jd(f„ ^/\g\e{{ip},t) = 0. (20) 

Geometrically speaking, this condition is fulfilled if the manifold is embedded in such 
a way that the center of mass of the surface parametrized by e coincides with the 
origin. For example, the sphere in Fig. lb is centered. 

The framework above will be illustrated by several examples in the subsequent 
chapters. According to our previous discussion, we follow the recipe as outlined before: 

(i) Parametrization of the manifold by a convenient set of generalized coordinates; 

(ii) Calculation of the Jacobian as well as the metric tensor; 

(iii) Derivation of the inverse metric tensor and noise amplitudes; 

(iv) Deduction of the Fokker-Planck and Langevin equation from the general expres¬ 
sions discussed in this section; 

(v) Solution of the model either analytically using Fokker-Planck equations or nu¬ 
merically via Langevin simulations. 


4 Isotropic self-propelled motion 

The director e(t) parametrizes points on a hypersphere S‘^~^ for self-propelled motion 
in d spatial dimensions. We choose |e(t)| = 1 for convenience. In this section, we start 
out with the case d = 2, derive the equations of motion in three spatial dimensions 
and then conclude with a general discussion for higher dimensional cases. 


4.1 Self-propelled motion in d=2 

In this subsection, we discuss the motion of a SPP in d = 2 which was studied 
in [32,39]. In two spatial dimensions, the director which determines the direction of 
motion of a particle describes diffusion on a circle of unit radius. Hence, the direction 
of motion can be parametrized by a single angle ip according to Eq. (4). The Langevin 
dynamics of ip is well-known, cf. Eq. (5). As an exercise, we will deduce this expression 
using our general approach discussed in the previous section. Since the director is 
described by a single variable, Greek indices can be omitted in this context. From 
Eq. (4), we obtain the Jacobian 

( d cos ip \ 

1 = 

d sin (p I 

~dP~ J 

and the corresponding metric tensor 

g = J^-J=l, 
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which is just a number in this case. Hence, we obtain expression (5) by inserting 
the metric tensor in the general expression (19). The corresponding Fokker-Planck 
equation is obtained from (15) by inserting the metric, or directly from (5): 

dtpi(p,t) = Ddlp{(p,t). ( 21 ) 

The solution of this Fokker-Planck equation is obtained from separation of variables 

p{ip,t\ipo,to) = 7 ^ + -Y] cos[m((,?- (22) 

ZTT TT 


where t > Iq. The correlation function of the director reads 

p2TT 

{e{t) ■ e{t')) = / d(p cos{(p — (fi') p{(p,t\(p',t') = ^ . (23) 

Jo 

From the Taylor-Kubo relation (3), we obtain the correlation time as well as the 
diffusion coefficient 


Tc ^ , 


2D 


(24) 


which were derived in [32, 39] first. The correlation time determines the crossover 
timescale at which the persistent ballistic motion of a particle turns into random 
diffusive motion. 


4.2 Self-propelled motion in d=3 


We follow the recipe described above in order to derive the equations of motion for a 
SPP in three dimensions which is less trivial than the previous example as shown in 
the following. In three spatial dimensions, the director is parametrized by two angles 
9 and p via 

( ei\ /sin0cos(p\ 

62 = sin0sin(p , (25) 

63 J \ cos 6 j 


where p G [0, 27r) and 9 G [0, tt) in our convention. At first, we calculate the Jacobian 


J = 


/3ei 

de 

dei \ 

dip \ / 

de2 


de 

des 

\ de 

IS ■ i 

dip) ^ 


cos 9 cos p 
cos 9 sin p 
— sin 9 


— sin 9 sin p\ 
sin 9 cos p 
0 ) 


as well as the metric tensor 


g=J^.J = 



(26) 


The metric tensor has a diagonal form, i.e. we can use (19b) in order to derive the 
Langevin dynamics of the angles: 


(S) 

(S) 


d9 

dt 

d(p 

dt 


Vwut) + 


D 

tan 9 ’ 


sin 9 


Ut)- 


(27a) 

(27b) 
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The corresponding Fokker-Planck equation reads 


dp{d,(p,t) 

dt 


_9 


Dp 
tan 9 


D 


d'^p D d'^p 
89“^ sin^ 9 dip 


,2 ■ 


(28) 


By integration over the polar angle ip, we obtain a Fokker-Planck equation for the 
marginal probability density p{9, t) 


dp{9, t) 
dt 

with its time-dependent solution 
sin 9 




Dp 
tan 9 


D 


d'^p 

W 


p{9,t\9o,to) = '^{2m+ l)Pm{cos9)P^{cos9o)e 


(29) 


(30) 


m—0 


We denote Legendre polynomials by Pm- Equation (30) is the probability density 
function to find a particle moving with a angle 9 versus the ea-axis given that it 
started moving into the direction 9q at t = to- Since the motion is isotropic and 
the dynamics is Markovian, we can always assume that the initial direction of motion 
equals the direction e = (0,0,1) without loss of generality and use Eq. (30) to compute 
the correlation function: 

(e(t) • e(t')) = f d9 cos9p{9,t\0,t') = \- (31) 

Jo 

Note, that the correlation function is exponentially decreasing as well, i.e. the behavior 
of a SPP in d = 3 is similar to the motion in d = 2. However, the correlation time in 
two dimensions is a factor of two larger than in three dimensions. Accordingly, the 
diffusion coefficient is reduced by a factor of three: 


Tc = 


1 

TD ’ 


= 


At 

6D 


(32) 


4.3 Self-propelled motion in arbitrary dimensions 

Surprisingly, the properties of SPPs in two and three spatial dimensions are qualita¬ 
tively very similar. Their characteristics like the diffusion coefficient and velocity cor¬ 
relation time differ by prefactors. Naturally, the question arises whether the behavior 
of SPPs is qualitatively similar in all spatial dimensions or whether a crossover/critical 
dimension exists above which the qualitative behavior changes. In this section, we an¬ 
swer this question by showing that the velocity correlation function of a SPP is indeed 
exponentially decreasing irrespective of the spatial dimensionality. Furthermore, we 
use this opportunity to present another approach for the description of the diffusion 
of the director. 

The parametrization of the director by angles as discussed in previous sections 
becomes increasingly difficult in higher dimensions. For completeness, we present the 
stochastic dynamics of the angles in appendix A. Here, we apply a different approach 
to study the dynamics of the director in d dimensions using Cartesian coordinates [40]: 


(I) 


dcj 

dt 


dcj 

dt 


d 

-D{d - l)ej -I- ^ {Sij - eiCj) ^j{t), 

d 

'i/‘2D ^ ] {5ij — CiCj) ^j{t)- 

i=i 


(33a) 


(S) 


(33b) 
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This parametrization is completely equivalent to the parametrization by angles which 
can be proven by explicitly carrying out the change of variables from {cj} to angles 
{(/?^}. However, care must be taken that the rules for the transformation of variables 
from ordinary calculus can only be applied if the stochastic differential equation is 
interpreted in the sense of Stratonovich [35], i.e. in Eq. (33b). In this context, we stress 
again the fact that both stochastic differential equations in (33) describe the same 
stochastic process. In particular, the length of the director |e(t)| = 1 is conserved in 
both cases^. They do only appear different at a first glance because their interpretation 
is different. Moreover, the dynamics described by (33) is Markovian. 

In the following, we will calculate the correlation function of the director in ar¬ 
bitrary dimensions. For this purpose, we multiply (33a) by e(f'), where t' < t, and 
average over many realization of the process. Therefore, we obtain a linear equation 
for the director correlation function 


d(e(t) ■ e{t')) 
dt 


-D{d-1) (e(t).e(O). 


This equation is readily solved with the appropriate initial condition. We obtain: 


(e(f)-e(t')) = e-^(‘^-i)|*-‘'l. 


Hence, the correlation function decreases exponentially in all spatial dimensions d. 
We obtain the correlation time by integration, cf. Eq. (7), 


1 

“ D{d-1) 


(34) 


in this general case and the diffusion coefficient of a SPP in d dimensions: 


V 


X 


Dd{d-1)' 


(35) 


Consistently, these general results reproduce our earlier results. Thus, we have shown 
that the qualitative properties of isotropic self-propelled motion do not depend on 
the spatial dimensionality. The only difference is the dependence of the correlation 
time on the spatial dimension. 


5 Anisotropic self-propelled motion 

We aim at generalizing the isotropic SPP model to account for anisotropic environ¬ 
ments. We stress that anisotropic diffusion is fundamentally different from a biased 
motion due to some external forcing which trivially induces biased particle flow and 
ballistic motion at large timescales rather than diffusion. Here, we focus on SPPs 
with an anisotropic directional persistence which may be due to an anisotropic sub¬ 
strate. For the sake of clarity, we restrict our discussion to the two-dimensional case. 
However, the generalization to d = 3 is straightforward from the general framework 
discussed in previous sections. 

Let us assume that the director moves on a manifold parametrized in polar coor¬ 
dinates in the following way^ 

e=M=r(^)f^°®^Y (36) 

\eyj ysmifj’ ysmifj’ '■ ’ 

^ In particular, the term Ci oc —D{d — l)ei in Eq. (33a) ensures the conservation of the 
length |e| = 1 if the corresponding equation is interpreted in Ito sense. 

® The choice of this parametrization a priori excludes a number of curves in two dimensions 
which are not relevant in the context of self-propelled motion. 
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Fig. 2. Visualization of the director parametrized by 
(36). Exemplarily, we have chosen an ellipse with ec¬ 
centricity e. Furthermore, we represent the arc length 
I by a blue line, cf. (40). 


where r(ip) is a 7r-periodic function'^. Curves parametrized in this way are distorted 
circles, e.g. an ellipse as illustrated in Fig. 2. The ellipse is parametrized by 


r{(p) 


1 

— COS^ (/? ’ 


where the parameter e G [0,1) denotes the eccentricity. The limiting case of a circle 
(isotropic motion) corresponds to e = 0 and r{ip) = 1. 

The angle ip determines the direction of motion of the SPP in two dimensions. 
The derivation of the dynamics of p is straightforward from the general discussion in 
section 3. We only need to calculate the metric tensor which reads: 


g{p) = [r{p)]^ 


dr{p) 

dp 


(37) 


The concrete form of the dependence r{p) is of minor importance, since the Langevin 
and Fokker-Planck dynamics for the angle p and the probability density p(p,t), re¬ 
spectively, depend on g{p) only. 

We derive the Langevin equation for the angle p from (19): 


(S) 


dip 

dt 



(38) 


We conclude that the anisotropy is reflected by a state dependent noise amplitude 
{anisotropic directional persistence). If the particle moves in a direction p where 
g{p) is large (preferred direction of motion), the level of fluctuations is reduced. In 
contrast, fluctuations are amplified in regions where g{p) is small. Note, however, 
that the particle moves at constant speed. The anisotropic motion arises from the 
non-homogeneous distribution of times a particle moves in a certain spatial direction. 

The corresponding Fokker-Planck equation for the dynamics of the probability 
density p{p,t) is obtained from the Langevin equation (38) for the angle: 


dp{^, t) ^ r^d_ 
dt dp 

Without further knowledge, this equation is difficult to solve. However, our geometric 
interpretation of the problem allows us to derive the time dependent solution of the 

If r{p) is not yr-periodic, the center of mass of the curve would be shifted away from zero. 
Consequently, this would imply biased ballistic motion. By assuming yr-periodicity, equation 
(20) is obeyed. 


d_ I P{P.t) \ 


\/9{‘p) V / 


(39) 
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Fokker-Planck equation without specification of g{(p). In order to do so, we exploit 
the fact that our problem can be mapped to the one-dimensional diffusion equation 
which can be easily solved analytically. In a first step, we parametrize the manifold 
in terms of the the arc length [43] 

1= [ dip' ^g{p '), (40) 

Jo 

cf. Fig. 2 for an visualization, giving this change of variables a transparent interpre¬ 
tation. We denote the total length of the curve by 

L= dp' V gi^p') ■ 

Jo 

We introduce the probability density p{l, t) to find a particle within the length interval 
\l,l + dlj. By substitution of variables in (39), we indeed find ~ as expected - the 
ordinary diffusion equation 


dp{l,t) ^ d'^p{l,t) 

dt 5/2 ■ 

The time dependent solution of the diffusion equation in one dimension with periodic 
boundary conditions, p{l + L,t) = p{l, t), is found by separation of variables. In fact, 
it has already been discussed in the context of isotropic self-propelled motion in two 
dimensions. Using this solution and subsequent resubstitution of p for the arc length 
yields the time dependent solution of Eq. (39): 


p{p,t\po,to) = 


Vgip) 


1 + 2J2 


COS 


m—1 


2'Km 


d<y?' '/gW) 


^—{2'KmjL)‘^D (t — to) 


(41) 


Notice that p{p,t\po,to) yf p{pojt\p,to) in contrast to what we learned for isotropic 
diffusion, c.f. (22). The stationary distribution po{p) = p{Pit\po,to —oo) is given 
by geometric properties only 


»(*=) = (42) 

since the strength of angular fluctuations D determines the timescale of relaxation, 
but it does not influence the stationary distribution as argued in the context of equa¬ 
tion (6). The distribution of the arc length po(0 is a constant, 

Po{l) = p (43) 

solely determined by the total arc length L. 

These results illustrate our general discussion in section 3. SPP motion, as it is 
considered in this work, is a combination of two processes: (i) motion at constant 
speed along the director and (ii) Brownian motion of an imaginary particle on a 
manifold described by the director. In turn, the imaginary particle is found with 
equal probability on every point of the manifold in the long-time limit since the 
dynamics of this particle is Brownian. This is reflected by expression (43). However, 
the distribution of the direction of motion p may be anisotropic as indicated by 
expression (42) depending on the way the manifold is embedded in space. The central 
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Fig. 3. Angular Distribution po{'^) for different values of the anisotropy a. Lines represent 
theoretical results, cf. (49), and points show results of numerical Langevin simulations of 
the particle dynamics. Parameters in numerical simulation: #particles N = 10 000, At = 
4 X 10-^ D = 1, no = 1. 


quantity which relates the motion of the imaginary particle on the manifold to the 
direction of motion of the SPP is the metric tensor g{ip). 

The general solution given by equation (41) allows us to estimate dynamic quan¬ 
tities as well, such as the characteristic correlation time. For large observation times 
t ^ to, the relaxation of the angular probability distribution is determined by the 
slowest mode: 




Vajp) 

L 


1 + 2 cos 


27r 




„-( 27 r/L)"D(t-to) 


(44) 


Thus, the characteristic relaxation time reads 


4,'k'^D 


(45) 


The perimeter of a circle equals L = 2tt confirming the prediction of isotropic self- 
propelled motion in two dimensions, cf. (23): = D~^. Note, however, that for a 

deformed circle whose perimeter is kept constant at L = 2 tt, the overall relaxation 
time does not change but is still determined by 

(46) 


Thus, the correlation function is expected to decay exponentially with a character¬ 
istic timescale Tc, when initial transients are neglected. The length of the transient 
depends on the timescale separation between the first and the second mode in (41). 
Interestingly, the ratio of characteristic times of the leading and the next-to leading 
order mode does not depend on model parameters but is a fixed value: 

—V/f— 

2t:VdJ / V47r/D 



This time scale separation is large enough to be observed experimentally or in simula¬ 
tions and justifies our reasoning above. The estimation of the correlation time allows 
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Fig. 4. Mean squared displacement as a function of time for two different anisotropies: 
a = 0.1 (left) and a = 0.9 (right). The full mean squared displacement as well as the 
components parallel and perpendicular to the preferred axis are shown. The dotted line 
indicates the limiting behavior <^(Zlr(t))^^ ~ 4:T>xt with the approximate diffusion coefficient 
'Dx = cf (47). Parameters in numerical simulation: ^particles N = 10 000, At = 

4 X 10■^ D = 1, no = 1. 


us to evaluate the diffusion coefficient via the general relation (8): 


V 


X — 



2D' 


(47) 


Interestingly, the diffusion coefficient does not depend on the degree of anisotropy. 

In order to visualize our results, we consider an illustrative example where the 
metric is given by 


9{<f) = 


1 2 n 2 

1 — a cos (f 
1 - 0/2 


(48) 


The degree of anisotropy is measured by the parameter a G [0,1]. For a = 0, the 
metric tensor is equal to a constant and isotropic motion is expected. If a is larger than 
zero, the metric tensor possesses maxima at = -k 12 + 1117: and minima for </_ = utt 
where n G Z. According to Eq. (38), angular fluctuations are suppressed if the particle 
moves along the j/-axis of the laboratory frame, whereas fluctuations are increased if 
the particle moves parallel to the x-axis. Thus, the diffusion of the SPP is anisotropic. 

We denote the mean squared displacement along the x-axis by ((Z\rj_(<))^ ) and along 

the y-axis by where the indices T and || indicate the mean squared 

displacement perpendicular and parallel with respect to the preferred direction of 
motion (y-axis). 

We have chosen this example (rather than the ellipse shown in Fig. 2) because the 
resulting equations are compact and therefore easier to interpret. In Fig. 3, we plot 
the stationary angular distribution 


Po{+) = 


1 — a cos^ ip 

(2 — a)7r 


(49) 


for different anisotropies. The dependence of the arc length parameter I on the polar 
angle reads in this case: 

a sin(2(p) 


l{ip) = ip 


2 — a 


2 


(50) 
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Fig. 5. Velocity correlation function {v(t) • v(0)) obtained from particle based Langevin 
simulations of the microscopic dynamics (38). The correlation function is averaged over all 
initial directions of motion. After a short initial transient, the velocity correlation decays 
exponentially with a characteristic time which is determined by the slowest mode in the 
probability density function (41): Tc = D~^. Parameters in numerical simulation: ^particles 
N = 500000, At = 4:X 10-^, D = 1, no = 1. 


Measurements of the mean squared displacement as well as the velocity correlation 
function from numerical Langevin simulation of (38) are depicted in Fig. 4 and Fig. 5, 
respectively. Interestingly, the long time behavior of the mean squared displacement 
does not change due to the anisotropy in accordance with (47). However, one can 

observe a splitting of the parallel (^|z\r|| and perpendicular component of 

the mean squared displacement with increasing anisotropies, see Fig. 4. Moreover, the 
predictions on the long-time behavior of the correlation function and the characteristic 
correlation time agree well with the prediction (46), cf. Fig. 5: the correlation function 
does possess exponential tails with a characteristic time determined by Tc = D~^. 


6 Summary &. Outlook 

We presented a geometric view onto SPP motion at constant speed. In order to do 
so, we divided the problem into two separate parts: (i) motion at constant speed in 
the direction given by the director in the d-dimensional physical space; (ii) stochas¬ 
tic dynamics of the director itself which takes place on a curved, abstract, {d — 1)- 
dimensional compact manifold. We saw that process (i) and (ii) are intimately linked. 
For example, the mean squared displacement can immediately be obtained from the 
correlations of the director. This implies that all relevant information is contained 
in the dynamics of the director which can be mapped to Brownian motion - char¬ 
acterized by a noise strength D - on a compact manifold. One essential step of this 
approach has been the parametrization of the manifold using generalized coordinates 
or angles. We showed that the dynamics of each generalized coordinate is described 
by a Langevin equation with a stochastic and a deterministic term, both proportional 
to D. We used this formalism to indicate that SPPs perform isotropic diffusion in the 
d-dimensional physical space if the manifold corresponds to a hypersphere retrieving 
previous results [32,39]. We also showed that the SPP performs anisotropic diffusion if 
the manifold does not correspond to a hypersphere, but to the surface of an arbitrary 
surface embedded in d-dimensional physical space such that Eq. (20) is obeyed. In 
summary, the proposed approach provides an unified framework to model isotropic 
as well as anisotropic diffusion of SPPs in any dimension. 
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A remark on the origin of anisotropic diffusion is in order here: The presence of 
a force field induces, asymptotically, ballistic motion, while we describe SPP with 
an anisotropic directional persistence remaining asymptotically diffusive. These two 
problems - anisotropic diffusion and random motion in a force field ~ are fundamen¬ 
tally different. Anisotropic diffusion is often caused by an anisotropic substrate. For 
example, the motion of eukaryotic cells moving on a pre-patterned surface exhibits 
a preferred axis of motion [33,34]. A similar reasoning may apply to gliding bacteria 
where, for instance, trails of polysaccharides on the substrate can induce anisotropic 
motion of the bacteria [44]. 

Finally, the proposed geometric approach to SPP motion provides a simple frame¬ 
work to describe the motion of cells, bacteria or other microorganisms in isotropic as 
well as anisotropic environments. It may be viewed as a coarse-grained description of 
the actual system under consideration in the sense that a complex self-propelled ob¬ 
ject moving on an anisotropic substrate is modeled by a point particle with anisotropic 
directional persistence. In this regard, it is essential that model parameters are im¬ 
mediately linked to observables which are easily accessible experimentally. We have 
illustrated this statement by a two dimensional example studied in section 5. The sta¬ 
tionary probability distribution to find a particle moving in a certain direction sheds 
light on the anisotropy of angular fluctuations. Similarly, the amplitude of these fluc¬ 
tuations enters into both the diffusion coefficient and the characteristic relaxation 
time of the velocity correlation function. Once these quantities have been measured, 
the point particle model can caricature the motion of a real-world active particle 
sufficiently well. 


A Diffusion on n-spheres 

For completeness, we state the equations of motion of a Brownian particle diffusing 
on a n-sphere. This surface, embedded in d = n -I- 1 dimensions, can be parametrized 
by n angles. We use the convention [45] 


Cl = cos ipi, 

n-i 

= cos(/?^ sin(/?m, = 2,d - I 

m—1 


d-1 

Cd = U S\TLLpra, 

m—1 

where V5i,...,d-2 G [0:'^) and S [0,27r). The metric of the sphere in this coordi¬ 
nates is diagonal and reads 


9p.<^ 


= S„ 


1 


/X = 1, 


\nl=\sinVm /x = 2,3,...,d-1. 
Thus, we obtain the following Langevin equations for the angles [41]: 
d(pi 


(S) 

(S) 

(S) 


^ =v^m + D^, 

at tan (fi 


dt 




mil sin A 

d(pd_i 


■ + D 


d — fj, — 1 


tan niii sin^ 


nliiSin:p,„ ^ 


^i = 2,3, ...,d-2, 
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This set of equations in combination with Eq. (1) are convenient to simulate self- 
propelled motion in d spatial dimensions. 
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